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Abstract 



One-loop quark contributions to the dimension-four gluon condensate term in the operator product 
expansion (OPE) of the scalar glueball correlation function are calculated in the MS scheme in the chiral 
limit of ri{ quark flavours. The presence of quark effects is shown not to alter the cancellation of infrared 
(IR) singularities in the gluon condensate OPE coefficients. The dimension-four gluonic condensate 
term represents the leading power corrections to the scalar glueball correlator and, therein, the one-loop 
logarithmic contributions provide the most important condensate contribution to those QCD sum-rules 
independent of the low-energy theorem (the subtracted sum-rules). 

The QCD correlation function of scalar gluonic currents 



is used to study the properties of scalar gluonium via QCD sum-rule techniques ^H]. The current J(x) 
is the lowest-order version of the operator fi(a)G 2 (x), which is renormalization-group invariant for chiral 
quarks H] . As first noted in E] , the one-loop gluon condensate contribution to (JJ) 



where {ao, bo, b\} are numerical coefficients (see (|26 j) -(|28 [l below), represents the leading condensate con- 
tribution to those sum-rules independent of the low-energy theorem .7 a , and so provide important non- 
perturbative effects within sum-rule analyses of scalar glueballs. 

The one-loop coefficients &o an d b\ were evaluated in [H] in the absence of quark effects {i.e. the rif = 
limit). As these rif = gluon condensate effects have been used in a number of sum-rule analyses where 
the effects of three quark flavours have been included in the perturbative part 0, it is necessary to extend 
the results of [5] to enable self-consistent sum-rule analyses in the presence of rif chiral quarks. 
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(1) 



J{x) = aG" (x)G a ^{x) = aG 2 {x), 



(2) 




(3) 



Since the gluonic current © is gauge invariant, the operator product expansion (OPE) of the corre- 
sponding two-point operator contains only those local operators which are gauge invariant, equations of 
motion, or BRS variations jlUj : hence, for massless quarks, the relevant OPE is given by 



i J d D xe iqx TG 2 (x)G 2 (0) 
= l(q 2 )+C(q 2 ) :G 2 (0): +V(q 2 )q°q X 



:GV(0)G"1(0):4to:G 2 (0): 



+ equation of motion & BRS invariant operators + higher dimensional operators (4) 

where I{q 2 ), C(q 2 ), T>(q 2 ) are the Wilson (OPE) coefficients and where the two colons indicate normal 
ordering. For notational convenience, the normal ordering symbol and the spacetime argument will subse- 
quently be omitted from the right-hand side of (|4"|). 

The scalar gluonic correlator Q is obtained from (j3J) through multiplication by a 2 followed by a vacuum 
expectation value (VEV). Physical matrix elements of the equation of motion and BRS invariant operators 
vanish, as does the VEV of the gluonic operator proportional to T>(q 2 ). Therefore, up to dimension four 
operators, the sole contributions to Q stem from I(q 2 ) (perturbation theory) and C(q 2 ). This article is 
concerned with the latter. 

As in , the nonzero momentum insertion technique (NZI method) [SJ |5] is employed to compute OPE 
coefficients. This method sandwiches the OPE between two external single gluon states with momenta k\ ^ 

permitting the simple separation of (non-physical) operators whose VEV is zero (i.e. BRS invariants 
and equations of motion) from the physical operator G 2 , thus simplifying the operator-mixing effects 
originating from the renormalization of composite operators [I]. The NZI method also facilitates the 
analysis of infrared aspects of the OPE, since the nonzero difference between the external gluon momenta 
s = (ki — /C2) 2 provides an infrared regulator. As shown in [S] for the nf = case, all infrared logarithms 
(i.e. log(s), log (k 2 ), and log (fef)) cancel in the calculation of the OPE coefficients. This then allows the 
use of on-shell external gluon states 

\ n ) = | e n; k n , c n ), £ n • k n = k n = (5) 

to sandwich the OPE, immediately eliminating all non-physical operators. The colour index c n associated 
with the gluon states © is included for completeness, but only leads to trivial overall colour factors when 
taking matrix elements of colour singlet objects such as As will be discussed below, the IR cancellation 
in the OPE coefficients is not altered by the inclusion of quarks, and hence these simplifications can be 
applied even in the presence of chiral quarks. 

Sandwiching (j3J) between on-shell gluon states yields 

i J d D xe iq - x {l\TG 2 (x)G 2 (0)\2) = C(g 2 )(l|G 2 |2) +V(q 2 )q^q\l\G a lia G a \ - ^ 9(7 xG 2 \2) 

+ contributions from higher dimension {d > 4) operators, (6) 

where the matrix elements of all non-physical operators have been eliminated through the use of on-shell 
external gluon states. To simulate the effects of a VEV, the direction averaging operator J dq which, for 
example, leads to the identity 



/ 



d<? q a qpf{q 2 ) = j^q 2 g a pf(q 2 ), (7) 



is applied to ©• Consequently, the term proportional to T>{q 2 ) is annihilated (as in the VEV). In addition, 
it should be noted that (1|G 2 |2) = 0(k 2 2 ) whereas the remaining terms on the right-hand side of © go like 
0(kf 2 )- in this way, contributions relevant to the computation of C(q 2 ) are easily identified. Therefore, 
Eq. © implies 

i J dqd D x e iq - x (l\TG 2 (x)G 2 {0)\2) = C(q 2 )(l\G 2 \2) + 0{k\ 2 ). (8) 
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Thus, at leading and next-to-leading order (respectively denoted by the (0) and (1) superscripts in what 
follows) in the bare fields, we have the set of equations 

i J dqd D x ^ x {l\TG 2 B {x)G 2 B {m) m = C(°)(g 2 )(l|G||2}(°) + G(kl 2 ), (9) 
i J dqd D x e iq x (l\TG 2 B (x)G 2 B {0)\2)W = (q 2 )(l\G 2 B \2)^ + C^\q 2 ){l\G 2 B \2)^ + 0{k\ 2 ), (10) 

where the subscript B denotes a perturbative expansion in terms of bare (unrenormalized) quantities. 
Consider first © which, in terms of Feynman diagrams, is given schematically by 



\ dq{ 



which implies 




+ crossed diagram > = C^\q 2 ) | 



} + 0{k\ 2 



(11) 
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[(Ail • k 2 )(ei • e 2 ) - (h • e 2 )(k 2 • ei)] + 0{ki <2 ) 

= 4C(°)(g 2 )5 ClC2 [(h ■ fcXei • e 2 ) - (fe • e 2 )(A :2 • e x )] + 0(kf >2 ). (12) 



Therefore 

C(°)(g 2 ) 



16 



(13) 



and we note that, at leading order, C(g 2 ) is unaffected by the inclusion of chiral quarks. The spacetime 
dimension D must be kept arbitrary until later stages of the calculations. 

Quark contributions to C(q 2 ) do, however, show up at next-to-leading order, but only through that 
contribution to the left-hand side of (fTU|) stemming from the diagram (and its crossed counterpart) depicted 
in Figure ^ This diagram does not generate any IR singularities in the OPE analysis since the external 
gluon momenta regulate the IR behaviour associated with the quark loop. Hence, the conclusion that 
all IR singularities cancel in the calculation of the gluon condensate OPE coefficient [5] is upheld in the 
presence of chiral quarks. As in 0, this cancellation then permits an expansion in the external gluon 
momenta prior to evaluating Feynman integrals (as opposed to an expansion after evaluating Feynman 
integrals used to explicitly show the IR-cancellation) . 




Figure 1: Feynman diagram containing quark effects contributing to the amplitude on the left-hand side 
of pi). 



Having justified an expansion in external momenta prior to evaluating Feynman integrals, it immedi- 



ately follows that 1 



(1|G||2)W 




= 0(k{ 2 ), (14) 
as a consequence of the resulting massless tadpole integrals, and hence (fTT)|) can be simplified to 

i Jdqd D xe^ x (l\TG 2 B (x)G 2 B (0)\2)^ = 4S c ^C^(q 2 )[(kx ■ k 2 )(e x ■ e 2 ) - (h ■ e 2 )(k 2 ■ ex)]+0(kl 2 ). (15) 

Calculation of the left-hand side of ()15|) constitutes a sizeable project. Fortunately, the vast majority of 
the required work has already been completed in Therein, all contributions within the framework 
of purely gluonic QCD were considered. Therefore, to extend the result to include chiral quarks, only 
those additional diagrams which admit internal quark loops need to be summed. As previously noted, at 
next-to-leading order, there are, in fact, only two such diagrams. Thus, 



i dqd D x e iq - x (l\TG 2 B (x)G 2 B (0)\2)^ 



mc 2 




- — } Pi • fcXei • 6 2 ) - (h ■ e 2 ){k 2 ■ e x )\ + O(fc^) 



+ crossed diagram > 
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+ 1 T " T nf 



x [(h • fc 2 )(ei • e 2 ) - (h ■ e 2 )(k 2 ■ ex)] + G{kU) 



(16) 

(17) 
(18) 



where D = 4 + 2e, 4 = - — 7 B + ln(47r) and where the first term on the right-hand side of (|16|) is the purely 
gluonic contribution calculated in [SJ. Together, Eqs. (|15|) and (|18|1 imply that 



- + log 

e 



Lastly, recalling (j2j) and using (|T3|) and (|19|) . we find 

i J d D x e iq x (0\TJ B (x)J B (0)\0) 

^ )(g 2 )+C«(g 2 )l (a B G%) + ... 

+ log 



+ 



131 13 
~6 



(19) 



I 4 a 

4 «b i p + 



B 
7T 



m 

6" 



B 

11\ t I 
1/ \ 



+ 



131 13 
"24" ~ 36 



(20) 



where the dots represent contributions to the scalar gluonic correlator arising from operators of dimension 
other than four. 



1 One-particle reducible self-energy contributions to the external gluons are ignored on both sides of COJ- 
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Eq. (|20j) is expressed in terms of bare quantities and so must be renormalized. Renormalization of 
composite operators is, of course, complicated by operator mixing [3]. Briefly, renormalized versions (no 
subscript) of a and G 2 are defined as 

a = Z~ 1 a B , (21) 
G 2 = Z G 2G 2 B + • • • (22) 

where the dots in (|22j) represent contributions from equation of motion and BRS invariant operators. Vac- 
uum expectation values of these omitted operators vanish and so they do not contribute to the calculation 
at hand; only Z G i is actually required. In the MS scheme the renormalization constant is |3] 

a ( 11 nt\ 1 , , 

Zo2=Z „ = 1 + _(_-_i)_ ; (23 ) 

thus, aG 2 is renormalization-group invariant at next-to-leading order. An analysis of the renormaliza- 
tion of G 2 using the NZI method is presented in jHJ, where Zq2 is easily distinguished from similarly 
defined renormalization constants corresponding to the omitted operators; hence Z G 2 can be determined 
by consideration of only a single amplitude: (1|G^|2). Noting that 

4 e , „ 

D = l -r < 24 > 

substitution of ((2"Tj) and 1)22(1 into (|2"U)l yields the following result in the MS scheme: 

i J d D xe iq - x (0\TJ(x)J(0)\0) 
= perturbation theory 



In {1 + - 

7T 



"i- l-M , {-<?\ , ^49 5 \]\ 2 (25) 

' {aG ) 



+ higher dimension condensate contributions. 
Finally, recalling @ allows us to identify the coefficients 

(26) 
(27) 

(28) 



«0 


= 4a, 


'49 


b 


= 4a | 






,12 


bi 


= 4a | 


' m 







appearing in Eq. (j2J). Modification of these one-loop results for a change of operator basis to f3{a)G 2 from 
aG 2 (in the operator J and/or in the OPE) can be achieved by algebraic rearrangements. 
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